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Abstract 

Aims: The hedgehog signaling pathway plays an important role in EMT of pancreatic cancer cells, but the precise 
mechanisms remain elusive. Because S100A4 as a key EMT moleculer marker was found to be upregulated upon Glil in 
pancreatic cancer cells, we focused on the relationship between Shh-Glil signals and SI 00 genes family. 

Methods: On the base of cDNA microarray data, we investigated regulating mechanism of Glil to some members of S100A 
genes family in pancreatic cancer cell lines firstly. Then, the regulation of Glil to S100A4 gene was studied by molecular 
biology assays and the pro-metastasis effection of Glil -dependent S100A4 was investigated in vitro. Finally, the expressions 
of Shh, Glil, S100A4 and E-cadherin in pancreatic cancer tissues were studied by using immunohistochemistry assays. 

Results: Five members of the S100 genes family, S100A2, S100A4, S100A6, S100A11, and S100A14 were found to be 
downregulated significantly upon Glil knockdown. Glil enhancer prediction combining with in vitro data demonstrated 
that Glil primarily regulates S100A family members via cis-acting elements. Indeed, the data indicate S100A4 and vimentin 
genes were upregulated significantly by Shh/Glil -expression increasing and E-cadherin was significantly reduced at the 
same time. Migration of PC cells was increased significantly in a dose-dependent manner of Glil expression (P<0.05) and 
siS100A4 significantly reversed the response of PC cells induced by L-Shh transduction (P<0.01). 

Conclusion: Our data establish a novel connection between Shh-Glil signaling and S100A4 regulation, which imply that 
S100A4 might be one of the key factors in EMT mediated by Shh-Glil signaling in pancreatic cancer. 
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Introduction 

Pancreatic cancer represents one of the leading causes of 
cancer-related mortality in industrialized countries [1]. The poor 
prognosis of this disease is mainly due to its early systemic 
metastasis [1-3]. A large number of studies have shown that the 
epithelial mesenchymal transition (EMT) might be a key 
mechanism of pancreatic cancer cells detaching from the primary 
tumor site and spreading to distant organs. However, the EMT- 
initiating signaling pathways remain elusive now in pancreatic 
cancer cells. Recently, hedgehog (Hh) signaling pathway has been 
demonstrated to be an important promoter of tumor growth in 
several gastrointestinal tract cancers [4,5]. As to pancreatic cancer, 
it has been demonstrated that aberrant activation of Hh signaling 
pathway was resulted from sonic hedgehog (Shh) overexpression in 
the majority of cases [6-8]. A large number of studies have shown 



that the main mechanisms of Hh signaling pathway in cancer cells 
were to promote epithelial mesenchymal transition (EMT) process 
[9-11]. Moreover, it has been demonstrated that blocking this 
pathway significantly inhibited metastasis of tumor cells in vivo 
and in xenograft models [12,13]. 

Since target genes of Glil were the key mechanisms of Hh 
signaling pathway, we attempted to understand its pro-metastatic 
mechanisms by identifing the targets of Gli 1 in pancreatic cancer 
cells. As a key pro-metastatic target gene of Glil in pancreatic 
cancer it must have the following characteristics: First, there are 
effective Glil cis-acting elements in its gene sequence; Second, its 
transcription was regulated positively by Glil; Third, it was 
upregulated in the majority of pancreatic cancer tissues and its 
expression level was positively correlated with Hh signaling; 
Fourth, its must be a key pro-metastatic functional factor. In 
previous study, we have performed a systematic research about the 
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Table 1. The primer sequences for real time RT-PCR assays. 



Gene 


primer sequences 


Glil 


F: 


5'-TCTGCCCCCATTGCCCACTTG-3' 




R: 


5'-TACATAGCCCCCAGCCCATACCTC-3' 


Shh 


F: 


5'-CGGAGCGAGGAAGGGAAAG-3' 




R: 


5'-TTGGGGATAAACTGCTTGTAGGC-3' 


Patchedl 


5'-CGGCGTTCTCAATGGGCTGGTTTT-3' 


5'-GTGGGGCTGCTGTTTCGGGTTCG-3' 


GAPDH 


F: 


5'-ACGGATTTGGTCGTATTGGG-3' 




R: 


5'-TGGAAGATGGTGATGGGATT-3' 


E-cadherin 


F: 


5'- CAATGCCGCCATCGCTTAC -3' 




R: 


5'- CAAAATGCCATCGTTGTTCACT -3' 


Vimentin 


F: 


5'- CAATGAGTCCCTGGAACGCC -3' 




R: 


5'- CACGAAGGTGACGAGCCATT -3' 


S100A2 


F: 


5'- CGGTCCAGGATGCCCAGTC -3' 




R: 


5'- GGCTCCCAGGGTGAGGATTTAT -3' 


S100A4 


F: 


5'- ACTCGGGCAAAGAGGGTGA -3' 




R: 


5'- CTGGGCTGCTTATCTGGGAA -3' 


S100A6 


F: 


5'- GGGGAGACTCGTCACCAGGC -3' 




R: 


c /— / — rf~t~ A At — iv — [— r/— /— at/— A /-/— / — r n' 
j - GG I LL AAG I L I I LL A I LAGLL I -3 


S100A10 


F: 


5'- GCCGCCTCCCTCTACCCAC -3' 




R: 


5'- GCCTTTATCCCCAGCGAATTT -3' 


S100A11 


F: 


5'- AGGGCGTGGGTTGAGGAGA -3' 




R: 


5'- TGTGGAGATGATGACAGAAAGGC -3' 


S100A13 


F: 


5'- GCTGGGAGGAGCGGTTAGA -3' 




R: 


5'- GCAAAGGTGAAGAAGGTGGTGA -3' 


S100A14 


F: 


5'- CAGTGTCGGTCAGCCAACG -3' 




R: 


5'- TCCCAGAAACTCCTGAACTCCA -3' 


S100A16 


F: 


5'- TGAACTGGGGTCCCTTTGTG -3' 




R: 


5'- TCAGGACGGACCCAGAATCA -3' 


S100P 


F: 


5'- AAGAGGCTGCCAGTGGGACA -3' 




R: 


5'- TTGGTCAGGGTCTGCGTGC -3' 


TCHH 


F: 


5'- TTCGGAGACCACATGACCCTA -3' 




R: 


5'- GCTCTTCCCGTTCTTGCCA-3' 



doi:1 0.1 371 /journal.pone.0096441 .t001 



target gene profiles upon Gli 1 in high-metastatic pancreatic cancer 
cell line through cDNA microarray and found that 5 members of 
this gene family were upregulated by Glil. Especially, the S100A4 
as a key moleculer marker promoting EMT process in pancreatic 
cancer was found to be upregulated more than 3 times in this 
study. On the base of these studies we focused on the relationship 
between Shh-Glil signals and SI 00 gene family. 

In this study, we analyzed the Glil biding sites within DNA 
sequences of SI 00 genes family with bioinformatics tools and 
databases. Moreover, we attempted to identify new connections of 
Shh-Glil signaling pathway with SI 00 genes family and to 
demonstrate pro-metastatic function of S 1 00A4 gene mediated by 
Shh-Glil signals in pancreatic cancer. 

Materials and Methods 

Cell cultures 

The human PC cell lines, BxPC3, AsPC-1 and Panc-1, were all 
commercial cell lines and we obtained these cell lines from the 



Institute of Cytology Chinese Academy of Science.These cell lines 
were widely used in pancreatic cancer research [14-18]. Three cell 
lines were all cultured in RPMI-1640 supplemented with 10% 
fetal calf serum at 37°C in a humidified atmosphere of 5% C02 in 
air. 

Vectors construction and cells infection 

Lentiviral transfer vectors for human Glil shRNA or Shh 
cDNA were constructed by Genechem Co., Ltd, Shanghai, China. 
This system includes the lentiviral vector pLVTHM, the envelop 
plasmid pMD2G, and the packaging plasmid pRsv-REV and 
pMDlg-pRRE. The lentivirus-Shh (L-Shh) contained a 3.3-kb 
coding sequence for Shh and the lentivirus-Glili (L-Glili) 
contained Glil small hairpin RNA (5'-CTCCACAGGCATA- 
CAGGAT-3') to the targeting sequence of the shRNA as 
previously described [19]. The lentivirus-control (L-C) which did 
not include the Glil interference sequences or Shh cDNA 
sequences served as control. Finally lentiviral constructs were 
verified via DNA sequencing to ensure accuracy. Recombinant 
lentivirus was produced by transient transfection of 293T cells 
following a standard protocol. When BxPC3, AsPC-1 and Panc-1 
cells were about 50% confluent in RPMI-1640 containing 2% 
FCS, they were infected with the lentivirus constructs at MOI 5. 
The cells were harvested after 72 hours for further experiments. 
To identify functional L-Shh and L-Glili constructs, we routinely 
analyzed Shh and Glil expression by qRT-PCR and western- 
blotting. 

S100A4 Transient Knockdown 

We used a RNAi sequence (siS100A4) and a negtive control 
sequence (siS 1 00A4 MIS) which has shown efficient knockdown of 
S100A4 in human colon cancer HT29 cells in a previous report 
[20]. S100A4 knockdown was monitored by determining its 
mRNA and protein levels after 48 h upon siRNA transfection, 
respectively. And then the cells were used to transwell assays. Cells 
were lipofected with the siRNAs using Lipofectamine 2000 
(Invitrogen) in according to the manufacturer's instructions. 

RNA extractions and real time RT-PCR assays 

Total RNA samples were extracted with Trizol reagent 
(Invitrogen) according to the manufacturer's protocol and lOOng 
total RNA was reverse transcribed in 20 |il volume and 2 |xl 
cDNA used for PCR according to the manufacturer's instructions. 
(TaKaRa Biotechnology, Dalian, China). The primer sequences 
were seen in table 1. CT (cycle threshold) values were standard- 
ized to CT values of GAPDH. 

Protein extractions and western-blotting assays 

Total protein samples were extracted with RIPA buffer 
according to the standard method and the samples were 
normalized for protein content by a commercially available kit 
(Bio-Rad Company). Equal amounts of protein were used for 
western-blotting assays according to the standard protocol for Shh, 
Glil, S100A4, E-cadherin, VIM (vimentin) and GAPDH detec- 
tion. 

Transwell assays 

Cell invasion assay-24 samples kits (Chemicon, Bedford, MA) 
were used to study the invasion/ migration of PC cell lines. After 
24h incubation, migration counts of PC cells were obtained by 
photographing the membrane through the microscope. Counts 
were taken in 5 highest cells areas at a high power magnification 
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Table 2. The primers for XChIP assays. 



Gene primer sequences 

S100A4 F: 5'-GCTGTGGCACTTACCGCATC-3' 

R: 5'-CCTCCGTCTTACGTGCATGTG-3' 
S100A6 F: 5'-GCCCTCCTTGCCCACATTG-3' 

R: 5'-TCTTCGTGACACGTGACTCGG-3' 
S100A2 F: 5'-GTGCTGAAGTGTTGACTGAAGGG-3' 

R: 5'-GGCAGGAGAATGGCGTGAA-3' 
GAPDH F: 5'-TACTAGCGGTTTTACGGGCG-3' 

F: 5'-TCGAACAGGAGGAGCAGAGAGCGA-3' 

doi:1 0.1 371 /journal.pone.0096441 .t002 

(x200). The mean value of the fields counted was considered as 
the migration count of PC cells. 

The Formaldehyde Cross-linking Chromatin 
Immunoprecipitation (XChIP) assays 

The XChIP assays were carried out as the previous studies [21]. 
In brief, the chromatin of AsPC-1 cells (3xl0 7 ) was collected with 
1 mL IP buffer containing protease inhibitor cocktails. Chromatin 
was sheared by using a sonicator in an ice box (6 rounds of 10s 
pulses, 350 W, 60 second intervals) and Crosslinking was reversed 
by adding 20 uL of 5 M NaCl overnight at 65°C. DNA was 
extracted using phenol/chloroform assay. 20 |iL of DNA was 
electrophoresed on a 1.5% agarose gel and the rest was used as 
INPUT DNA. ChlP-grade goat polyclonal Glil antibody (0. 1 ug/ 
mL) (Santa Cruz Company) was used for immunoprecipitation, 
the mouse IgG (0. 1 ug/ mL) (Santa Cruz Company) was added as 
a random control, RNA polymerase II antibody (0.1 |J.g/mL) as a 
positive control and [i-actin antibody (0.1 |ig/mL) as a negative 
control. The primer of the promoter region of GAPDH (Santa 
Cruz Company) was used as a negative control. (Table 2). 

Luciferase reporter vectors assays 

The S100A4 promoter luciferase reporter vectors were 
constructed by Genechem Co., Ltd, Shanghai, China. Briefly, a 
1.5-kb cDNA fragment (from -766 to +734 of human S100A4 
gene) with Xhol and Hindlll sites was synthesized and cloned into 
the Xho I and Hind III sites of the pGL3 basic luciferase vector to 
produce pGL3-1.5 S100A4 containing three Glil biding sites 

[ACCCACCAC (-349 357), TGGGTGGTG (-227 235) 

and CTGGTGGGG (126- 134)]. In the mutagenesis vectors, Glil 
potential binding sites sequences were replaced respectively by 
GATTCTTAA sequence that has no known binding sites for any 
transcription factors [22]. The single Glil biding site mutagenesis 
vectors included pGL3-1.5 S100A4 Mut 1 [GATTCTTAA (- 

349 357), TGGGTGGTG (-227 235) and 

CTGGTGGGG (126-134)], pGL3-1.5 S100A4 Mut 2 [AC- 
CCACCAC (-349 357), GATTCTTAA (-227 235) and 

CTGGTGGGG (126-134)] and pGL3-1.5 S100A4 Mut 3 

[ACCCACCAC (-349 357), TGGGTGGTG (-227 235) 

and GATTCTTAA (126-134)]. The multiple Glil binding sites 
mutagenesis vectors included pGL3 1.5 S100A4 Mut 1-2 

[GATTCTTAA (-349 357), GATTCTTAA (-227 235) 

and CTGGTGGGG (126-134)], pGL3 1.5 S100A4 Mut 1-3 

[GATTCTTAA (-349 357), TGGGTGGTG (-227 235) 

and GATTCTTAA (126-134)] and pGL3 1.5 S100A4 Mut 2-3 

[ACCCACCAC (-349 357), GATTCTTAA (-227 235) 

and GATTCTTAA (126-134)]. The final vector constructs were 



verified via DNA sequencing to ensure accuracy. AsPC-1 cells 
were transfected with 5 |Xg of either the pGL3-1500 S100A4 or the 
pGL3-1500 S100A4 mutant constructs and 2 ug of the renilla 
luciferase construct (pGL4.75, Promega) as an internal control for 
normalization using the Lipofectin reagent method protocol. 
48 hours after transfection, AsPC-1 cells were used for the 
luciferase assays according to the protocol from the dual luciferase 
reporter kit (Promega) and the volues were readed using a 
luminometer. Each experiment was repeated at least 3 times, each 
time in triplicate. 

Patient samples 

In this work, 102 banked paraffin-embedded cancer tissues from 
different PC patients (from 63 men and 39 women, mean age 56.7 
years, range 44—75 years) were analyzed by immunohistochem- 
istry. All patients have agreed this research. All samples were 
obtained from the tenth people's hospital of Tongji university, 
China and all patients have signed a consent document, in which 
they agreed that the surgical excision of malignancy would be used 
in medical research before their operations. The ethics committees 
of the Tongji Universities approved the consent procedure and the 
studies. 

Immunohistochemistry assays 

The paraffin-embedded PC tissues were used for identification 
of Shh, Glil, S100A4 and E-cadherin. Shh polyclonal antibody 
was purchased from R&D Company and the others (Glil, S100A4 
and E-cadherin) were from Sant Cruz Company. The immuno- 
histochemistry assays were carried out as the preceding studies and 
the expression levels of Shh, Gli 1 , S 1 00A4 and E-cadherin genes 
were graded based on previously described guidelines [23]. 

Statistical Analysis 

For all statistical analyses, a software package SPSS 17.0 (SPSS, 
Inc, Chicago, IL, USA) was used. Continuous variables were 
expressed as means ± SE and a non-paired Student's t-test was 
used for statistical evaluation. The correlations between expression 
of Shh, Glil, S100A4 and E-cadherin proteins and the relation- 
ship between the clinical characteristics of the patients and the five 
proteins were analyzed using the Spearman rank test. P<0.05 was 
considered to be statistically significant. 

Results 

Differential expression of S100 gene family upon Glil in 
AsPC-1 cells 

To delineate the gene expression induced by Shh-Gli 1 signals in 
pancreatic cancer cells, AsPC-1 cells were used to cDNA 
microarray assays comparing lentivirus-control vs lentivirus-Gli 1 i 
cells. When focusing on the EMT-related genes in cDNA data, we 
found the positive relationship between Shh signals and SI 00 
genes family[21]. In our screen (the expression signals included 23 
members of SI 00 genes family except for S100A17 and S100A18) 
10 genes, S100A10, S100A11, S100A13, S100A14, S100A16, 
S100A2, S100A4, S100A6, SI OOP and TCHH, were expressed 
and 13 genes, S100A1, S100A12, S100A3, S100A5, S100A7, 
S100A15, S100A8, S100A9, S100B, S100G, S100Z, RPTN and 
FLG were not expressed in AsPC-1 cells. Comparison analysis 
from lentivirus-control group vs lentivirus-Glili group cells showed 
that S100A2, S100A4, S100A6, S100A11 and S100A14 were 
upregulated upon Glil. (Figure 1A). The ten genes, S100A10, 
S100A11, S100A13, S100A14, S100A16, S100A2, S100A4, 
S100A6, SI OOP and TCHH, were picked for validation by real- 
time RT-PCR. As showed in Figure IB, S100A2, S100A4, 
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Representative Public ID Gene Symbol Fold Change 



A 



K 1 n it A AT A7D 

NM_005978 


O A A A A O 

S100A2 


2.548 


NM_002961 


O A A A A A 

S100A4 


3422 


NM_0 14624 


O A A A A d 

S100A6 


2.839 


t I r ,1 A AAAA A 

NM_002966 


O A A A A A A 

S100A1 0 


1.325 


NM_005620 


A*( AA A A 

S100A1 1 


1.898 


NM_005979 


S100A13 


1.217 


NM_020672 


S100A14 


2.096 


AA0451 84 


S100A16 


1.175 


N 30878 


TCHH 


1.031 


NM_005980 


S100P 


0.654 



B 




□ L-GM1i 
■ L-SHh 



J.^i.ri .ri.ri 



S100A2 S100A4 S100A6 S100A14 S100A11 S100A10 S100A13 S100A16 TCHH S100P 

Figure 1. Expression of the SI 00 gene family. A: cDNA microarray data about S100 gene family; B: The differential expression levels of partial 
members of SI 00 gene family by qRT-PCR. The expression of GAPDH is as a control. * P<0.05, ** P<0.01. 
doi:1 0.1 371 /journal.pone.0096441 .g001 



S100A6, S100A11 and S100A14 were found to be upregulated, 
with no significant changes in the other five genes. These data are 
consistent with that obtained from the cDNA microarray. 

Predicted Gli1 biding sites and putative Gli1 target genes 
within S100A gene family. 

To further study the mechanism of Shh-Gli 1 signaling pathway 
on the regulation of the SI 00 gene family, we downloaded the 
DNA sequences of S 1 00A genes from Gene Bank and analyzed 
homologous sequences of Glil biding site (GACCACCGA) by 
utilizing Blastz bioinformatics tool. The bioinformatics data 
showed that almost all members of S100A gene family contain 
many highly conserved homologous sequences (the difference is no 
more than one base) and these putative Glil biding sites exhibited 
the characteristics of the cluster arrangement. The Table 3 
showed the putative Glil biding sites within the proximal 
regulatory region containing intergenic regions, intragenic regions, 
untranscripted regions and untranslated regions (no more than 
5.0 kb away from TSS). The data probably confered a somewhat 
higher overall degree of evidence that these candidates are in fact 



differentially expressed in the setting of pancreatic cancer, possibly 
as downstream targets of an aberrantly reactivated. 

The new target gene and its effective Glil binding sites 
identification in AsPC-1 

The PCR primers for XChIP assays were designed according to 
1.5kb-long promoter sequences of S100A2, A4 and A6 genes, 
which containing several predicted Glil biding sites (64~72, 

326-334, 418-426 for S100A2; -349 357, -227 235, 

126-134 for S100A4; -246 254, 786-794 for S100A6. the 

homology of each site is 89%.). (Table 3). The results of the 
XChIP assays showed that transcription factor Glil bound to the 
promoters of S100A2, 4 and 6 genes respectively in the cultured 
AsPC-1 cells (Figure 2). 

Results of luciferase assays showed that promoter activities of 
pGL3-1500 S100A4 increased with expression levels of Glil gene 
in a dose dependent manner (P<0.05). The pGL3-1.5 S100A4 
Mut 2 or 3 had no significant effect on luciferase activity (P>0.05), 
however, pGL3-1.5 S100A4 Mut 1 decreased luciferase activity 
significantly compared with pGL3-1.5 S100A4 (P<0.05), which 
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Table 3. The putative Glil biding sites in S100A gene 


family. 






Gene 


Site 


Sequence 


Homology 


S100A1 


350-358 


ATGGTGGGT 


89% 




1415-1423 


CTGGTGGGT* 


100% 




1780-1788 


GGGGTGGTC 


89% 


S100A2 


-4932 4940 


TGGGTGGTA 


89% 




-3703 3711 


ACCCCCCAG 


89% 




-3323 3331 


TGGGTGGGC 


89% 




-2870 2878 


AGCCACCAG 


89% 




-2366 2374 


CTGGGGGGT 


89% 




64-72 


TGGGTGGGC 


89% 




326-334 


CTGGTGAGT 


89% 




418-426 


CTGGTGGGA 


89% 


S100A3 


-4900 4908 


ACCCTCCAG 


89% 




-3205 3213 


GACCACCCA* 


1 00% 




-1567 1575 


ACCCACAAG 


89% 




491 -499 


TGGGTGGTG 


89% 


S100A4 


-349 357 


ACCCACCAC 


89% 




-227 235 


TGGGTGGTG 


89% 




126-134 


CTGGTGGGG 


89% 


S100A5 


-947 955 


GTGGTGGGT 


89% 




424-432 


CTGGTGGGT* 


100% 


S100A6 


-246 254 


CTGGTGGGG 


89% 




786-794 


TGGGTGGTT 


89% 


S100A7 


-4794 4802 


ATGGTGGGT 


89% 




-4276 4284 


TGGGTGTTC 


89% 




-4141 4149 


CACCACCCA 


89% 




-3686 3694 


CACCACCCA 


89% 




-2717 2725 


CTGGTGGCT 


89% 




-99 107 


ACCCACCTG 


89% 


S100A8 


-4865 4873 


ACCCACCAC 


89% 




-3944 3952 


ACCCACCTG 


89% 




-218 226 


GACCACCAA 


89% 




484-492 


TCCCACCAG 


89% 


S100A9 


-4428 4436 


GACCACCCC 


89% 




-3295 3303 


TGGGTGGGC 


89% 




-3169 3177 


ACCCACCTG 


89% 




-2294 2302 


CTGGTGGGG 


89% 




-793 801 


CTGGTGGGC 


89% 


S100A10 


-3501 3509 


ACCCACTAG 


89% 




383-391 


ACCCACCCG 


89% 




522-530 


TGGGTAGTC 


89% 




617-625 


CTGGTGGGG 


89% 




784-792 


TTGGTGGGT 


89% 




2638-2646 


TACCACCCA 


89% 








89% 


S100A11 


-2685 2693 


CGAACACCCA 


89% 


S100A12 


-2510 2518 


TGGGTGCTC 


89% 




-2048 2056 


CTTGTGGGT 


89% 




-225 233 


CGGGTGGGT 


89% 




78-86 


TAGGTGGTC 


89% 


S100A13 


- 1 340 1 348 


AGGGTGGTC 


89% 
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Table 3. Cont. 





Gene 


Site 


Sequence 


Homology 




-1305 1313 


ACCCACCTG 


89% 




3400-3408 


GGGGTGGTC 


89% 




3731-3739 


CTGGTGGGT 


100% 




4660-4668 


ATGGTGGGT 


89% 


S100A14 


-1656 1664 


TGGGTGCTC 


89% 




138-146 


TGGGTGGTT 


89% 


S100A15 


-1620 1628 


GTGGTGGGT 


89% 




1233-1241 


ACCCCCCAG 


89% 


S100A16 


-483 491 


CTGATGGGT 


89% 




110-118 


TGGGAGGTC 


89% 




270-278 


CTGGGGGGT 


89% 




448-456 


GGGGTGGTC 


89% 




2194-2202 


CTGGTGGGG 


89% 




2682-2690 


CACCACCCA 


89% 



doi:1 0.1 371 /journal.pone.0096441 .t003 



suggest that site 1 might be enhancers of Glil. The luciferase 
activity of plasmids contains site 1 increased significantly (P<0.01) 
and that of plasmids contain site 2 or 3 did not change significandy 
(P>0.05). (Figure 3A, B, C). The sequence of site 1, ACCCAC- 
CAC, had also been demonstrated to be recognized and sis- 
activated by Glil proteins previously [24,25]. 

Regulation of Shh-Gli 1 signals to expression of S100A4 
and to invasion/migration of PC cells through mediating 
S100A4 in vitro 

To further study the pro-metastasis function of Glil -derived 
S 1 00 A4 in pancreatic cancer cells, five groups of PC cells, L-Gli 1 i, 
L-C, L-Shh, L-Shh+siS100A4 MIS and L-Shh + siS100A4, were 
used to analize the regulation of Shh-Gli 1 signals to S100A4, E- 
cadherin and VIM genes by real time RT-PCR and western- 
blotting assays. And then the same five group cells were used to 



evaluate invasion/migration of PC cells in vitro by transwell 
assays. The results of qRT-PCR and western-blotting showed that 
expression levels of S100A4 and VIM genes were increased 
significandy by Shh/Glil -expression increasing. In contrast, the 
expression levels of E-cadherin were significantly reduced at the 
same time. (Figure 4A, C). Moreover, the S100A4 knocked-down 
signifiandy reversed downregulated E-cadherin and upregulated 
VIM induced by L-Shh transduction. (Figure 4B, D). The results 
of transwell assays showed the migration of PC cells were 
significandy decreased in L-Glili group and increased in L-Shh 
group respectively compared with L-C group (P<0.05). (Fig- 
ure 4E). Moreover, siS100A4 significantly reversed the response of 
PC cells induced by L-Shh transduction (P<0.01). (Figure 4F). 
Thus, it seems that S 1 00 A4 mediated by abnomal activated Shh- 
Gli 1 signaling pathway could be one of key pro-migration factors 
in pancreatic cancer cells. 
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Figure 2. Binding of Glil to promoter region of S100A2, 4 and 6 genes in AsPC-1 cells analyzed by XChlP-PCR assay. M: DNA Marker; 
PC: RNA polymerase II antibody for XChIP positive control; IP: Glil XChlP; IGG: mouse IgG for XChIP random control; NC: (3-actin antibody for XChlP 
negative control. 
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Figure 3. Identification of effective Glil binding sites on the S100A4 promoter. A: S100A4 luciferase activity increases with increasing 
expression level of SHh in AsPC-1 cells. The pGL3-1.5 S100A4 and Renilla luciferase vectors were transiently transfected into L-Gli1i infected, L-C 
infected or L-Shh infected AsPC-1 cells respectively. Results were normalized for transfection efficiency using Renilla luciferase and the L-C infected 
group was arbitrarily given a value of 1. B: Relative luciferase activity of the S100A4 promoter Glil binding site mutants. AsPC-1 cells transfected by L- 
Shh were transiently transfected 5mg of each reporter construct including pGL3-1.5 S100A4, pGL3-1.5 S100A4 Mutl, Mut2 and Mut3 respectively. The 
pGL3-1.5 S100A4 was arbitrarily given a value of 1 and the activities of the other transfections were adjusted relative to this activity. C: Site 1 was 
responsible for Glil transcription. Relative luciferase activity from different AsPC-1 cell groups, L-GN1i infection, L-Shh infection or L-C infection, 
transfected by different constructs, pGL3 1.5 S100A4 Mut 1-2, Mut 2-3 and Mut 1-3 with Renilla luciferase vectors. Results were normalized for 
transfection efficiency using Renilla luciferase and the L-C infected group was arbitrarily given a value of 1. * P<0.05, ** P<0.01. 
doi:1 0.1 371 /journal.pone.0096441 .g003 



Correlation between the expression of Shh, Glil, S100A4 
and E-cadherin proteins in PC tissues 

The results of immunohistochemistry showed that the positive 
stainings of Shh. Glil and S100A4 proteins localized almost 
exclusively to the neoplastic cells beside weakly positive staining of 
Shh in the islet cells and the chronic pancreatitic duct complex 
occasionally being seen in pancreatic tissues apart from cancer. In 
contrast, the protein expression of E-cadherin was markedly 
downregulated in PC tissues compared with normal tumor- 
adjacent tissues. The staining patterns of Shh and S100A4 were 
diffuse cytoplasmic type, that of E-cadherin was diffuse cytomem- 
brane type and that of Glil display perinuclear cytoplasmic and 
nuclear staining. Because the Glil was activated only after 
entering the nucleus, the samples with only cytoplasmic staining 
were counted as negative expression. (Figure 5). Immunoreactivity 
ratios were 78.43% for Shh (80 of 102), 51.96% for Glil (53 of 
102), 81.37% for S100A4 (83 of 102) and 48.04% for E-cadherin 
(49 of 102) in PC samples. The Spearman rank test analysis 
showed that significant positive correlations were noted between 
Shh, Glil and S100A4 at the protein level in PC (Shh vs Glil: 



CC = 0.697, P<0.01; Shh vs S100A4: CC = 0.476, P<0.01; Glil 
vs S100A4: CC = 0.510, P<0.01). Conversely, the significant 
negative correlations between the three proteins and E-cadherin 
were confirmed (Shh vs E-cadherin: CC = -0.278, P<0.01; Glil 
vs E-cadherin: CC = -0.207, P<0.05; S100A4 vs E-cadherin: 
CC = -0.285, P<0.01). 

The correlation between the expression of four proteins 
and clinicopathological features of the PC patients 

The statistical analysis data showed that expression levels of Shh 
(CC = 0.245, P<0.05 and CC = 0.254, P<0.05) and S100A4 
(CC = 0.265, P<0.01 and CC = 0.220, P<0.05) proteins were 
positively correlated with tumor size and lymph node metastases 
and that of Glil protein was positively correlated with lymph node 
metastases (CC = -0.370, P<0.01), while that of E-cadherin 
protein was negtive correlated with lymph node metastases 
(CC = —0.203, P<0.05). However, no significant correlation was 
found between the four proteins and age, gender, tumor location 
and differentiation. 
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Figure 4. Hh signaling pathway promoting invasion/migration of pancreatic cancer cells through mediating S100A4. A: The relative 
transcription levels of Shh, GN1, S100A4, E-cadherin and vimentin were regulated by L-Glil i/L-Shh transduction. B: The relative transcription levels of 
S100A4, E-cadherin and vimentin in L-Shh transfected cells were reversed by siS100A4 transduction. C: The expression levels of Gli1, S100A4, E- 
cadherin and vimentin proteins regulated by L-Glil i/L-Shh transduction. D: The expression levels of S100A4, E-cadherin and vimentin proteins were 
reversed by siS100A4 transduction. E: The invasion/migration of pancreatic cancer cells regulated by L-Glil i/L-Shh transduction were analyzed by 
transwell assays. F: The invasion/migration of L-Shh transfected cells were reversed by siS100A4 transduction. * P<0.05, ** P<0.01. 
doi:10.1371/journal.pone.0096441.g004 
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Figure 5. Expression of Shh, Glil, S100A4 and E-cadherin in 
pancreatic cancer tissues analyzed by immunohistochemistry. 

A: Strong cytoplasmic staining for Shh in pancreatic cancer cells (x400); 
B: Weakly positive staining for Shh in ductal complex of chronic 
pancreatitis tissues (x400); C: Weakly positive staining for Shh in islet 
cells of normal cancer side tissue (x200); D: Negtive staining for Shh in 
normal ductal epithelial/acinar cells of cancer side tissue (x200); E: 
Cancer cells show strong nuclear staining for GN1 and normal cancer 
side tissues were negtive (x400); F: Strong positive perinuclear 
cytoplasmic and part nuclear staining for Glil in pancreatic cancer 
cells (x400); G: Strong positive cytoplasmic staining for S100A4 in 
pancreatic cancer cells (x400); H: Positive cytoplasmic staining for E- 
cadherin in pancreatic cancer cells. Black arrows point to positive 
staining area and white arrows point to negtive staining area. 
doi:1 0.1 371 /journal.pone.0096441 .g005 

Discussion 

During tissue morphogenesis activation of Hh signaling 
pathway almost accompany the occurrence of EMT and this 
process is required for migration of Hh-responsive cells [26,27]. It 
was reported that the mechanisms of Shh-Gli 1 signaling pathway 
promoting malignant tumors involved in many aspects of 
malignant biological behaviors but its main role in tumors is to 
promote EMT and maintain tumor stem cells [28] . In this study, 
we provide a systematic research about the mechanism of Glil in a 
high-metastasis pancreatic cancer cell line, AsPC-1. The AsPC-1 
cell line contains both Ki-ras gene mutation and DPC4 gene 
deletion, which were typical pancreatic cancer molecular patho- 
logical change [29,30]. Moreover, the expression level of Glil was 
knocked down moderately by using the RNA interference system, 
which was more consistent with the actual molecular pathological 
change of pancreatic cancer than artificial gene knockout. So the 
target gene spectrum obtained in our screening was more 
consistent with the actual situation of pancreatic cancer cells. 



To date, at least 25 members have been reported to belong to 
SI 00 family in humans. Twenty-one of them are coded by genes 
clustered at chromosome locus lq21, known as the epidermal 
differentiation complex (EDC) involving in epithelial-derived cell 
differentiation [31]. It has been reported some members of this 
family were overexpressed in various cancer types and their 
functions might involved metastasis and EMT [32 — 34]. The EDC 
contains two groups of consecutively distributed S100A genes. 
One group contains S100A1, S100A13, S100A14, S100A16, 
S100A2-S100A7, S100A15, S100A8, S100A12 and S100A9 and 
the other contains FLG, HRNR (S100A18), RPTN, TCHH, 
TCHH1 (S100A17), S100A11 and S100A10. The other genes 
belonging to the subfamilies of S100B, SI OOP, S100Z and S100G 
are, respectively, located at chromosome loci 21q22, 4pl6, 5ql4 
and Xp22. It was been reported that S100A genes had a common 
origin in the molecular evolution while S100B, of SI OOP, S100Z 
and S100G had different origins, so S100A genes might have a 
universal conserved regulatory sequences [35]. Given their highly 
conserved Glil binding homologous sequences, it is speculated 
that they play similar role related to Gill, though we cannot 
exclude the possibility that they have different functions unrelated 
to Gill. Our results suggested that the expression of S100A genes 
family in PC cells has three manners: First, expressing and 
depending on Glil; second, expressing but not depending on Glil; 
third, not expressing. Combinating Glil enhancer prediction with 
microarray data, we speculated that Glil regulated SI 00 A genes 
through the cis-acting elements might be a fundamental mode of 
regulation. Previous study had reported that S100A7 and S100A9 
transcription were induced by GUI treatment in epidermal cells 
[25]. During evolutionary process some S100A genes expression 
might be turned off due to obtain unknown repressor sequence 
and then some of these genes might be turned on again but no 
longer were regulated by Glil as a result of obtaining additional 
cis-acting elements. This hypothesis may pardy interpret why the 
expression of SI 00 A family is characterized by tissue selectivity 
though it is still lack of sufficient evidence. Detailed comparative 
genomics studies between distinct species, as well as different 
SI 00 A genes might help to farther understand the exact 
mechanisms. In our microassay and the ensuing real-time PGR 
studies, several S100A genes, including S100A2, S100A4, S100A6, 
S100A11, and S100A14, were upregulated in response to Gill. 
Several of these S100A genes such as S100A2, A4 and A6 have 
been shown to be involved in pancreatic cancer. In this study the 
data of XChIP assays showed that transcription factor Glil bound 
to the regulatory sequences of S100A2, A4 and A6 genes 
respectively in AsPC-1 cells, which suggested that these genes 
might be cis-regulated by Glil. 

Studies have shown that the functions of SI 00 gene family 
members were complex and diverse and So far, only S 1 00A4 was 
found to be closely associated with EMT process. Moreover, in 
previous studies we have proposed the hypothesis that S100A4 
gene might be one of the key factors in EMT molecular network 
regulated by Shh-Glil signaling pathway in pancreatic cancer cells 
[19]. Therefore, S100A4 gene was used as an example to explore 
the relationship between Shh-Glil signals and SI 00 gene family. 
As a matter of fact, numerous studies have indicated that the 
S100A4 protein levels in cancer tissue, pancreatic juice and serum 
of patients with pancreatic cancer were significantly increased [36— 
38]. Some studies have implied that S100A4 could be used as a 
potential biomarker of pancreatic cancer [39]. However, why 
S100A4 gene was overexpressed selectively in PC was uncertain. 
Maybe S100A4 gene could be regulated by certain selectively 
activated signaling pathway at the same time in pancreatic cancer. 
In the present study, we for the first time found that the expression 
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levels of SI 00 A4 protein were correlated positively with activity of 
Shh-Glil signals in pancreatic cancer tissues and the Glil protein 
upregulated S100A4 mRNA through cis-activation manner in PC 
cells which established an exact molecular pathway from Shh-Gli 1 
signals to S100A4 in PC cells. 

Although activation of Hh signals almost accompanied the 
occurrence of EMT, to date, there was no evidence that Glil 
directly regulated Snail/ Slug transcription. It has been reported 
that S100A4 and E-cadherin were inversely regulated in several 
cell systems and S100A4 promoted the expression of the essential 
transcription factors, Twist and Snail, in the EMT process, as well 
as mesenchymal markers, including VIM and MMPs [40-42]. 
Our data demonstrated that the abnormal activated Hh signaling 
pathway promoted EMT at least partially through mediating 
S100A4 gene in PC cells. Since S100A4 was an important pro- 
metastatic factor in PC, the new connection of Shh with S100A4 
might be one of the key links within the molecular network during 
EMT of PC cells. 

In conclusion, these data suggested a model in which Shh-Glil 
signaling pathway promotes metastasis in pancreatic cancer by 
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